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Abstract— A generic semi-analytical approach is proposed for the self-consistent analysis of combinatorial 
frequency generation in stacks of binary nonlinear layers illuminated by a pair of pump waves. It is shown that the 
quasi-periodic stacks can be treated as the defected periodic structures with the defect located at the specific 
positions determined by a particular layer sequence. The developed technique combining the Harmonic Balance 
and Transfer Matrix Methods is illustrated by the cases of periodic and quasi-periodic (Fibonacci and Thue-Morse 
type) stacks of nonlinear dielectric layers. 
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Electric displacement in nonlinear anisotropic dielectric layers: 
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Quasi-Periodic Stacks 








Fibonacci stack of order g Thue-Morse stack of order q 
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+ Two TM plane waves of frequencies o, and o, are 
incident on the stacks at angles ©, and ©, 


* Nonlinearity is assumed weak and the three-wave 
mixing process is dominant 


+ The problem is linearized by harmonic balance 
method to obtain combinatorial frequency fields 
* [n the non-depleting wave approximation the fields 
are determined recursively at each frequency 
» Pump waves: only reflected and refracted fields 


» Combinatorial frequencies: excited & scattered fields 
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In the three-wave mixing process, fields H,(@) of TM waves 
satisfy Helmholtz equation ın each layer A and B at pump o» 


and combinatorial @,=@,+@, frequencies: 
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In a layer of type j=A,B in n™ primitive cell: 
+ At pump wave frequencies ©, 
(n) = n+ oi Z E oe -ic tik x WW 
H (o,)-| are +B" le , p-12 


* At combinatorial frequency (04 
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The problem has been reduced to evaluating 
amplitudes of waves refracted into a layer of 
type j = A, B in an n" primitive cell of the stack 
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where RK(Q,) - the reflection coefficient of the stack; 


n^. (o, ) - the transfer matrices of a subset of (n - 1) primitive 


cells preceding the layer of type j in the n™ primitive cell. 
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Amplitudes of the combinatorial frequency emission 
in the reverse (F,) and forward (F) directions 
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where My, (c SE a transfer matrix of the whole stack containing N, 
primitive cells 
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* The transfer matrices of the stack and any its 
subset depend on 
» stack configuration and layer sequence 
> Transfer matrices of the constituent layers 
+ The transfer matrices are to be calculated at each 
frequency: both pump waves and mixing products 


« The transfer matrices at frequency @, are used in 
the expression of F,, and $7 (o) 
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4 Transfer matrices for 


> Stack of n unit cell: ge = fp, = (Ma Loi lm, Loi II 
> 1“ layer in n unit cell stack: pe =ñ, M, (0) 
> Layers A and B: T. 0, ) and m,, Loi, 


+ Transfer matrices ol" have closed form relating 


i| With Bloch phase in periodic stacks 


* There 1s no advantage of using the closed form 
here — fields must be calculated inside each layer 


* Ihe recursive relation for 7, 7 7, , DL» Loi, Im, Loi ) | 
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Transfer matrices for Fibonacci & Thue-Morse 
stacks are defined by the recursive relations 


Fibonacci stacks of order 422 | Thue-Morse stacks of order q71 
$,719,4 J 5,2]: S,={A f, S,={AB} O,=19,.; Z Qi & Q 7-1Q Z Ou 


Transfer matrix Qo = {A}, Qo = {8} 
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MO. |= 40 5,10 m E xm 
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Fibonacci & Thue-Morse stacks can be 
treated like periodic stacks with defects 


Fibonacci stacks Thue-Morse stacks 
S; = {AB A'B AB A} Q; = {AB BA BA AB) 
Ss = {AB A'B AB A'B A'B} Q, = {AB BA BA AB BA AB AB BA} 
* Two types of primitive cells: * Two types of primitive cells: 
- regular (AB] - regular {AB} 
- "defective" {A'B i | - “defective” {BA} 
* Layer A in "defective" cells 1s a * Layers A and B in “defective” cells 
doublet A'ZAA with d", =2d A are interchanged 
* Both "regular and "defective" cells * Both "regular and "defective" cells 
have the layers in the same order have the same thickness d=d,+d, 


Positions of “defective” cells in the stacks are to be determined 
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* [he number of primitive cells 1n a stack of order q: 
Nr 
* The positions of the regular and defective cells are 
determined by their serial number n using the 
following recurrence relation at 7123 
EP 
t = i 
n (ray n odd 
2 


n— even 


n—|? 


where ¢, = 0 for the regular cell {AB} 
t, = 1 for defective cells {BA} 
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The number of primitive cells in a stack of order g: 
Poi! g 
N= ` 
®,,,-T,-1 
2 


,  q- even 
, q— odd 


where 6 - d, , * , » 1s Fibonacci number, 6, - &, =1 
The number of "defective" cells I: I, 2 0 at qx3 and at q24 
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Positions of additional A layers in Fibonacci stack 
of order q are defined by a row-matrix P, of length 
o, with I's in the columns for the first A layer of 
the doublets: 5S, = {AB AAB AB AAB AAB} 

P. =[00 100 00 100 100] 
e is defined by the recurrence relations: 

BB Loris Di = 

where u (o Si E 15.0, | Is a row-matrix with l's in P, column only; 
(a (o d - Q T ] is a Square Toeplitz matrix with |’s only at the 
secondary diagonal offset for ®, from the main diagonal. 
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Fibonacci 5; Periodic g=16 
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* The semi-analytical approach to modelling 
combinatorial frequency generation by periodic 
and quasi-periodic multilayers has been developed 


+ The technique provides a unified framework for 
the analysis of periodic and quasi-periodic stacks 


* It is shown that Fibonacci and Thue-Morse stacks 
can be treated like periodic structures with defects 


* The developed theory is illustrated by simulations 
and provides an insight in the mechanisms of 
three-wave frequency mixing in the binary stacks. 
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